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ABSTRACT 

Context. Analysis of thermal dust emission spectra for dust mineralogy and physical grain properties depends on comparison spectra, 
which are either laboratory-measured infrared extinction spectra or calculated extinction cross sections based on certain grain models. 
Often, the agreement between these two kinds of spectra, if available, is not yet satisfactory because of the strong influence of the 
grain morphology on the spectra. 

Aims. We investigate the ability of the statistical light-scattering model with a distribution of form factors (DFF) to reproduce 
measured infrared dust extinction spectra for particles that are small compared to the wavelength, i.e. in the size range of 1 yum and 
smaller. 

Methods. We take advantage of new experimental spectra measured for free particles dispersed in air with accompanying information 
O !■ on the grain morphology. For the calculations, we used DFFs that were derived for aggregates of spherical grains, as well as for 

compact grain shapes corresponding to Gaussian random spheres. In addition we used a fitting algorithm to obtain the best-fit DFFs 
for the various laboratory samples. In this way we can independently derive information on the shape of the grains from their infrared 
spectra. 

Results. With the DFF model, we achieve an adequate fit of the experimental IR spectra. The differences in the IR band profiles 
between the spectra of particulates with different grain shapes are simply reflected by different DFFs. Irregular particle shapes require 
a DFF similar to that of a Gaussian Random Sphere with cr=0.3, whereas roundish grain shapes are best fitted with that of a fractal 
aggregate of D/=2.4-1.8. The fitted DFFs generally reproduce the measured spectral shapes quite well. For anisotropic materials, 
different DFFs are needed for the different crystallographic axes. The implications of this finding are discussed. 
Conclusions. The use of this model could be a step forward toward more realistic comparison data in infrared spectral analysis of 
thermal dust emission spectra, provided that these spectra are dominated by emission from submicron grains. 

Key words, infrared: general — methods: data analysis — methods: laboratory — circumstellar matter — planetary systems: proto- 
C*"} , planetary disks 
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■ 1 . Introduction avoids the influence of an embedding material on experimental 

^7"; ' infrared extinction spectra. 
^ The analysis of mid-infrared dust emission spectra from stellar 
• *"H outflows, circumstellar disks, and other objects provides infor- These spectroscopic measurements of dust particles in 

^ ' mation about the dust mineralogy, grain sizes, and temperatures, aerosol also allow the investigation of the actual morphology of 

[ hence, about physical and chemical conditions in the respective the aerosol particles by filtering and subsequent scanning elec- 

environments. The information about the dust grain properties is tron microscope (SEM) imaging. Thus, the influence of morpho- 

based on the lattice vibrational bands of the dust particles that logical particle character istics o n the dust spectra c an be st udied 

mainly occur in the 10-50 fim wavelength range and that dom- in detail. Tamanai et al. d2006bl > and Tamanai et al. d2009bl have 

inate the thermal emission of warm dust. Detailed analyses of, reported a strong dependence of the measured profiles of the in- 

e.g., SpitzerIRS, and ground-based mid-infrared spectra includ- frared bands on grain shape and agglomeration. Although the 

ing interferometric data (VLTI/MIDI) have been published in the influence of the grain morphology on dust vibrational bands has 

last years for dust emission from accretion disks (e.g. van Boekel been known for many years (Bohren & Huffman I1983K these 

et al. l2004l l , debr is disks (e.g. Lisse et al. 120071) . and comets (e.g. effects have not been systematically studied in experiments be- 

Lisse et al. [2006t . fore. 

The analyses are often performed by applying^ 2 fits of a Tin- Furth ermore, Tamanai et al. d2006B demonstrated that the 

ear combination of either calculated or measured spectra of dust band fiJes dicted b the currentl used models often dif . 

grains to the observed spectra. Both approaches have recently fe considerabl from the experimentally measured ones. The 

been improved by developing (1) new theoretical models (Mm reasQn ^ ^ mQSt of thefje modelf , such afj Mig ^ CDE 

et al.|2005| Mm et aL l200^ tor calculating dust spectra and (2) a (Bohren & Huffman 1983 K and DH S (Min et al. I2QQ51 are not 

new experimental method for measuring spectra of particles dis- flexiWe gn h tQ tafce morphological effects sumcie ntly into 

persed in air (Tamanai et al.g006aj Tamanai et al.E006b), which account since M& h important for i nter p r etating astronomical 

spectra, we aim at (a) a clarification of typical properties of mea- 

Send offprint requests to: sured infrared band profiles depending primarily on the grain 
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Fig. 1. Synthetic particle shapes and DFFs used in the modeling for a) aggregates of spherical grains with different fractal dimensions 
D f and b) Gaussian Random Spheres with different standard deviations cr. 



shape, and (b) the investigation of the ability of an alternative 
theoretical model to take the shape effects into account. 

For this purpose, we compare a selection of aerosol- 
measured spectra of oxides and silicates for two different classes 
of grain morphologies with calculated spectra using the theo- 
retical approach of a form factor distribution (Min et al. 2006 ). 
Both the experiments used here and the theory are restricted to 
particles small compared to the wavelength. For infrared wave- 
lengths, this means a restriction to particles in the size range of 
approximately 1 yum and smaller. Consequently, the applicabil- 
ity of our results will be limited to interpreting emission (or ab- 
sorption) spectra from cosmic environments where such grains 
are present. If larger grains provide a strong contribution to the 
emission cross section, the spectrum will be influenced by grain 
size effects that cannot be treated by the present approach. This 
is especially the case when growth or radiation effects have re- 
moved submicron grains, such as in evolved protoplanetary disks 
and debris disks around luminous stars. 

The theoretical model and the particulate materials are intro- 
duced in the next two sections, before we compare the experi- 
mental and the simulated spectra in Sect. 4. Section 5 summa- 
rizes the results. 



2. DFF model for spectra simulation 

The distribution of form factors model (DFF model, Min et al. 
2006 ) is a statistical approach valid for ensembles of particles 
that are small compared to the wavelength (Rayleigh limit). It 
describes the averaged extinction (absorption, emission) cross 
section <C exl > per unit volume of such a particulate by the inte- 
gral over a continuous distribution of dipole polarizabilities de- 
pending on the dielectric function of the particulate material rel- 
ative to its environment e = s p /s e and a depolarization factor or 
form factor L, which is defined in the range 0<L<1 (Eq. 1). P(L) 
is the distribution of form factors over this range, which contains 



the information on the geometrical properties of the particles. 

V A Jo l/(e-l) + L 

Min et al. (2006) have demonstrated that a distribution P(L) 
(hereafter called the DFF) can be calculated for each particle 
shape, which gives a correct representation of the extinction 
cross section of that particle, provided that a spatial discretiza- 
tion with sufficient accuracy is possible. For a distribution of 
shapes, the individual DFFs have to be averaged. 

Another nice property of the DFF approach, which allows 
very intuitive understanding of the simulated spectra, is that ma- 
jor contributions to the integral occur when 

e ~ 1 - l/L (2) 

(so-called surface modes, see Bohren & Huffman [1983]). This is 
indeed the case in strong lattice vibration bands, where the real 
part of s takes a range of negative values (depending on wave- 
length), while the imaginary part is comparably small. Thus, re- 
lation (2) links the DFF directly to the extinction at correspond- 
ing spectral positions; e.g., for spherical grains (P(L)=<S(L= 1/3)), 
the extinction cross section spectrum will show a single res- 
onance at the wavelength where s=-2. Assuming a simple 
Lorentzian behavior of the dielectric function, the DFF at low 
L values is related to the strength of absorption at wavelengths 
longward of the sphere resonance (Re(e)<-2), whereas the DFF 
at larger L is related to absorption at shorter wavelengths (lower 
negative Re(e) values). We illustrate this in Fig.|2j and for a more 
precise treatment in the complex s plane see Min et al. (2006). 

The often used continuous distribution of ellipsoids model 
(CDE, Bohren & Huffman 1983) can be considered a special 
case of the DFF model with P(L)=2(1-L). Min et al. d2"0"06i ) have 
calculated DFFs resulting from certain particle shapes for the 
cases of (a) aggregates of spherical grains with different fractal 
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dimensions D / and of (b) particles in the form of Gaussian ran- 
dom spheres (GRS, Muinonen et al. [1996 ) with different stan- 
dard deviations <x of the surface modulation. We may note here 
that GRS grain shapes are used in scattering models for large 
grains as well (e.g. Volten et al. 2001). However, this is an en- 
tirely different (ray optics) approach, where the size distribution 
is also reflected by the (often extreme) GRS shape. In contrast, in 
our small-particle-limit approach the GRS model describes only 
the grain shape. 

The DFFs and illustrations of the particle morphologies are 
shown in Fig.Q] In both cases, the DFFs broaden strongly with 
variation of the respective parameter, i.e. when going from com- 
pact to fluffy (low-dimension) aggregates and from sphere-like 
to extremely structured GRS. 

However, the details are different in the two series. While 
the "main peak" of the DFF shifts to higher L with decreasing 
D f in the case of the spherical grain aggregates, the opposite is 
the case for the GRS with increasing <x. For the extreme cases, 
namely the aggregate with Dy-1.2 and the GRS with <x=0.7, 
these main peaks are located at L=0.6 and L=0.05, respectively. 
Consequently, we can expect the two types of DFFs to pro- 
duce band profiles peaking on opposite sides of the sphere res- 
onance position (Re(s) — 2). Figure [2] demonstrates this for less 
extreme geometries, i.e. an aggregate with D /=2.4 and a GRS 
with <r=0.3. The short-wavelength peak of the aggregate par- 
ticle spectrum (Re(e) — 1.5) corresponds to L~0.4, whereas the 
peak of the GRS spectrum (Re(e) — 6) corresponds to L~0.15. 
According to the DFF, it was expected to occur at L~0.1. The 
difference probably stems from the influence of Im(e), which in- 
creases towards longer wavelengths. 

For the aggregates of spherical grains, there is actually an- 
other peak in the DFF, which becomes stronger and stronger 
with lowering D /, namely at L=0. This component of the DFF 
produces a long-wavelength shoulder in the band profile in addi- 
tion to the short-wavelength peak (see Fig. The L=0 compo- 
nent indicates two- or one-dimensional elongated, i.e. plate- or 
needle-like, structures. Agglomerates of spheres with point-like 
contacts are in principle not expected to produce such a compo- 
nent. Its presence is a consequence of the limited resolution of 
the discretization, which does not allow between point-like con- 
tacts and extended contact areas characteristic of e.g. coalesced 
particle aggregates to be distinguished (cf. Andersen et al. 2006). 
This is a direct effect of the method used to compute the DFFs 
(based on the discrete dipole approximation). However, coales- 
cence of particles is often observed in real particulates, and its 
presence in the model distributions is not at all a shortcoming for 
their application, as we show. 

In this paper, we make use of these DFFs by comparing 
simulated infrared extinction spectra for such synthetic parti- 
cle shapes with the measured spectra of real particulates with 
known (different) shapes. We consider particulates composed of 
spinel, corundum, and forsterite (see next section), for which the 
wavelength-dependent optical constants s p (A) are known from 
the literature. We use data by Fabian et al. (2001a) for synthetic 
stoichiometric spinel, Querry et al. (1985) for corundum, and 
Sogawa et al. (2006) for forsterite. For the last two, the cross 
sections obtained for the individual crystal axes are averaged. 
For the environment, the dielectric constant s e is unity in the 
simulation of aerosol spectra, 2.31 in cases when the particles 
are embedded into KBr matrices, and 3.03 for the case of em- 
bedding into a Csl matrix. 
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Fig. 2. Extinction spectra of fractal aggregates with Dj — 2.4 
(solid black line) and Gaussian random spheres with <x = 0.3 
(dashed line) calculated for a Lorentzian oscillator type dielec- 
tric function. Re(e) is also shown (gray line). The vertical lines 
indicate the wavelength range where Re(e)=l-1/L can be ful- 
filled and the position of the sphere resonance (Re(e)=-2). 

3. Experimental spectra 

The experimental setup for the measurement of infrared extinc- 
tion spectra of dust particles dispersed in air (aerosol measure- 
ments) was described by Tamanai et al. (2006a 2006b, 2009a, 
2009b}. For comparison with the simulations based on the DFF 
model, we have selected two samples of each of the three 
minerals spinel (MgA^O/t), corundum (AI2O3), and forsterite 
(Mg2SiC»4), from these papers. The sample properties are listed 
inTableQ] 

The two samples of each material have significantly dif- 
ferent grain shapes. In general, one of the particulates can be 
characterized by a roundish shape or at least by having round 
edges, which is often the case for particles that have been con- 
densed from the gas or liquid phase in a relatively fast process, 
whereas the other particulate is characterized by an irregular 
grain shape with sharp edges, which is typical of particulates 
produced by crushing larger pieces of material, but may also re- 
sult from slow crystal growth. The spectra of the two classes of 
particulates differ significantly, as discussed by Tamanai et al. 
(2006b 2009b). Correspondence with astronomically observed 
infrared dust emission bands has been found to be better in some 
cases for the roundish shapes, especially for the oxides in AGB 
star outflows (Tamanai et al. |2009bl see also Posch et al. 119991 
and Fabian et al. 200la), but in other cases better for the irreg- 
ular shapes, especially for forsterite (Tamanai et al. 2006b[ see 
also Molster et al. 120021 Fabian et al. 2001b). Particle shapes and 
spectra are shown in Fig. [3] 

4. Results 

4.1. Simulated spectra vs. measured aerosol spectra 

Figure|4] shows the comparison of the simulated spectra for 
spinel, corundum, and forsterite particles with the measured 
spectra. The left column compares the spectra calculated for the 
aggregates of spheres to the measured spectra of the roundish 
particulates. The simulated band profiles consist in most cases 
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Fig. 3. Aerosol infrared extinction spectra and SEM micrographs of the spinel (left), corundum (middle), and forsterite (right col- 
umn) particulates. The solid lines are the spectra of the particulates with roundish shapes, which are shown in the images in the 
second row. The dashed lines correspond to the spectra of the sharp-edged grains (images in bottom row). The dark round filter holes 
in the image background are about 0.5 /mi in diameter except for the image of the roundish spinel sample (upper left micrograph), 
where they are 0. 1 fim in diameter. The numbers indicate SEM magnification. 

Table 1. Properties of the samples. 



Material 


Chemical formula 


Product info 


Processing 


Grain size 


Grain shape 


Reference 


Spinel 


MgAl 2 4 


Aldrich 




< 0.1 fim 


round 


T09b, CSp2 






Alfa Aesar 


Sedimentation 


< 1 fim 


irregular 


T09b, CSP1 


Corundum 


A1 2 3 


Alfa Aesar 




< 0.3 fim 


ellipsoidal 


T09b, CAC1 






GC Jena 


Sedimentation 


< 2 fim 


irregular 


T09b, CAC2 


Forsterite 


Mg 2 Si0 4 


Marusu 




< 0.3 fim 


ellipsoidal 


T06b, CF2 






Alfa Aesar 


Sedimentation 


< 1 /jm 


irregular 


T06b, CF1 



of peaks at short wavelengths and shoulders at longer wave- 
lengths (compare Fig. [2]). The latter strengthen with lowering 
fractal dimension of the aggregate, as do the L~0 components 
in the DFFs. The "primary" peak at shorter wavelengths, which 
weakens only between D/=2.8 and 2.4, remains present while 
shifting to shorter wavelengths. This often leads to a "rectangu- 
lar" band profile for low fractal dimensions, such as at D/<1.8, 
which actually reproduces the measured spectrum of the corun- 
dum sample well. The comparison with the other simulated spec- 
tra demonstrates clearly that in the case of these particles a 
strong long- wavelength band component is definitely required to 
fit the measured profiles at about 17 fim and 22.5 fim. The find- 
ing that a fractal dimension of Dy -1.8 provides the best match 
corresponds nicely to the observed agglomeration state of the 
corundum powder, which is mainly chain-like (see Fig. [3]). For 
the other samples of roundish grain shapes, the long-wavelength 
component is weaker but also present. For both the spinel and 
the forsterite samples, the D/=2.4 spectrum fits quite well to the 
peaks in the measured spectra. 



In the case of spinel, the region in between the two bands and 
longward of 20 fim shows an enhanced extinction that is not re- 
produced by any of the calculations. We interpret this feature as 
being caused by extinction from clumps of particles. Apparently, 
due to the small size of the individual grains, the clumps ap- 
pear compact enough to act partially similar to large grains. 
According to Mie calculations, grains of about 5 fim size already 
provide sufficient extinction (absorption and scattering) at these 
wavelengths. 

Moreover, we note that a few bands, such as the forsterite 
10 fim and 27.5 fim and the corundum 20 fim bands, are not 
satisfactorily reproduced by the same models as the majority of 
the bands. This is discussed in Sect. |4.2| 

The right column compares the spectra of the irregular par- 
ticulates with the simulations for the GRS of different surface 
modulation. As expected, the band profiles in these simulated 
spectra are dominated by peaks that shift continuously to larger 
wavelengths with an increasing standard deviation parameter. 
A shoulder remains at the short-wavelength edge of the pro- 
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files. The comparison with the measured spectra reveals excel- 
lent agreement for the cr=0.3 simulation in the cases of spinel 
and corundum (again apart from the /l=20yum band). For the 
forsterite, cr-0.3 is still the best-fitting calculated spectrum, even 
if it seems that a little lower cr may have given a better match. 
This is especially true for the lOyimi band, which we discuss sep- 
arately in Sect. 14. 21 

4.2. Fitting of DFFs and the problem of anisotropy 

An alternative way to compare theory and measurement is to fit 
DFFs to the measured spectra, hoping for a better reproduction 
of the measurements and aiming at the comparison of the ob- 
tained DFFs. We have developed an automatic fitting procedure 
that is able to handle anisotropic materials, i.e. contributions of 
up to three dielectric functions (in the case of forsterite, two for 
corundum) and even with the option of using different DFFs for 
each of these. In Sect. 14. 31 we show the result of this fitting for 
the irregular spinel. The fitting worked without problems for this 
isotropic material. 

In the previous section, we noted that a few bands in the 
spectra of the anisotropic materials show poor reproduction by 
the models, such as the 10 fim band of forsterite and the 20 
/im band of corundum. For reproducing them, a DFF peaking 
strongly between L=0.2 and 0.25 (close to that of GRS 0.1) 
would work best, but it would not reproduce the other bands. 
These bands belong to vibrational modes along certain crys- 
tal axes, so introducing different DFFs for different crystal axes 
could solve these problems. 

In general, deviations in the DFFs along different crystal axes 
are not an unlikely option, because the grains could be primarily 
elongated along certain crystal axes. Such grain shapes can be 
produced both by growth and shattering processes. Allowing for 
different DFFs for the different crystal axes lead to a satisfactory 
reproduction of the problematic bands both for forsterite and for 
corundum. Figure[5] shows as an example the fits for the irreg- 
ularly shaped forsterite with (a) only one DFF for all axes and 
(b) a separate DFF for the contributions of the vibrations along 
the crystallographic x-axis. The better match in the 10/im band 
is clearly visible. 

However, the DFFs for the different axes needed to be very 
different. The DFF required by the problematic bands were al- 
ways similar, with a narrow peak in the L=0.2-0.3 range. This 
changed only slightly for roundish grains, even if the main DFF 
peaked then on much higher L values. In Fig. [6] the two fitted 
DFFs are shown for the irregular forsterite particles (see Fig. [5]), 
where the difference is still quite moderate. The GRS model with 
<x=0.3, however, represents the case well for the y and z-axis 
DFF for the irregular corundum particles. Here, the difference 
from the x-axis DFF is very strong. Especially for the roundish 
grains, such extreme differences between the crystal axes and 
the nondependence on the grain-shape seemed unlikely. In some 
bands, we noted that the poor fit introduced by the additional 
DFF had to be compensated for by the other components. 

This pointed to a different reason for the mismatch in some 
bands, which turned out to be an incorrectness in the model, 
namely in the treatment of the anisotropy of the materials. The 
simple averaging of the absorption cross sections calculated for 
the crystallographic axes is correct for spheres that are small 
compared to the wavelengths and for ellipsoids, if the ellipsoid 
principal axes are associated with the crystal axes. When assum- 
ing complicated shapes, such as the GRS or aggregates, however, 
it is clear that this condition is not maintained. Here, the influ- 




Wavelength (|J.m) 



Fig. 5. Comparison of fits to the measured spectrum (gray line) 
of irregular forsterite particles using a single DFF for the contri- 
butions of all crystal axes (dashed line) and allowing two DFFs, 
i.e. a separate DFF for the contribution of the vibrational modes 
along the crystallographic x-axis (solid line). The DFFs obtained 
by the two-DFF fitting are plotted in Fig. [6] For comparison, the 
simulation using the GRS model with cr=0.3 is also shown (dot- 
ted line). 



Fit forst. irreg., y,z-axes - 
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Fig. 6. Comparison of the fitted DFFs for irregular forsterite par- 
ticles and irregular spinel particles (dashed line) with the GRS 
model with rx=0.3 (dash-dotted), which provides an excellent fit 
for the spectrum of the irregular corundum particles. 

ence of the polarization along one axis on the excitation along 
other axes has to be considered. 

This can be demonstrated by means of the discrete dipole 
approximation (DDA) model. The DDA allows one to take the 
full tensor nature of the refractive index into account. We will 
demonstrate the effects of this in a separate paper (Min et al. in 
prep). 

4.3. Reproduction of KBr and Csl spectra 

An interesting question is whether the DFF models can predict 
the influence of embedding media on the spectra of particulates. 
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Fig. 7. Comparison of measured (solid lines) compared to sim- 
ulated extinction spectra for roundish spinel particles embedded 
in different media. The DFF models used are the aggregates of 
spherical grains with D/=2.8 (dotted lines) and D/=2.4 (dashed 
lines). 



Tamanai et al. (2006b, 2009b) demonstrated the strong influence 
of potassium bromide and cesium iodide matrices used in the 
conventional pellet technique onto IR band profiles of embedded 
particulates. Within the DFF model, the embedding medium is 
taken into account simply via the relative dielectric constant e = 

Figure|7]and[8]show simulated spectra for spinel using iden- 
tical DFFs and different dielectric constants of the embedding 
medium {e e - \ for aerosol, 2.31 for KBr, and 3.03 for Csl) com- 
pared to corresponding measured spectra. For the irregularly 
shaped particles (Fig. 13, the DFF for GRS with cr=0.3 and a 
fitted DFF were used. The comparison shows that the former, 
which reproduces the aerosol measurements satisfactorily, gives 
a still quite imprecise representation of the measured spectra for 
embedded particles. However, this does not mean that the model 
fails to treat the influence of the embedding, because some ad- 
justment of the DFF (see Fig.|6]l leads to a nearly perfect fit of 
all the spectra measured in different matrices. 

Unfortunately, this is not the case for the roundish parti- 
cles (Fig. [7] see especially the 18/vm band). Here, the embed- 
ded particles are best matched by a DFF strongly dominated by 
the spherical character of the single grains or characterized by 
very compact aggregation, whereas for the aerosol spectra we 
had to use a DFF calculated for a lower fractal dimension of 
the aggregates. If this indeed reflects different morphologies of 
the spinel particles in aerosol and in embedded states, it could 
possibly indicate a compaction of the aggregates by the pellet 
pressing, or, since we have already noted that these spectra are 
also influenced by size effects, an aggregate size reduction by 
the thorough mixing with the matrix powder. The size of the ag- 
gregates of course could always be influenced by the preparation 
technique. Unfortunately, an inspection of the morphology of the 
embedded aggregates is not possible and a treatment of size ef- 
fects is beyond the capabilities of this theory. For the other sam- 
ples, which have larger grains and therefore show less clumping, 
we have not observed similar effects. 



Fig. 8. Same as Fig. [7] for irregular spinel particles. The DFF 
models used here are the Gaussian random spheres with cr=0.3 
(dashed lines) and a fitted DFF (see Fig.|6]l. 

5. Conclusions 

By comparison with experimental spectra, we have demon- 
strated the usefulness of the DFF approach for simulating the 
infrared extinction spectra of submicron-sized particulates com- 
posed of oxide and silicate minerals. Predefined DFFs with cer- 
tain characteristics can be used to achieve a satisfactory predic- 
tion for the shape dependence of the measured infrared band 
profiles. Good agreement between experimental and calculated 
band profiles was achieved by using 

1 . DFFs calculated for Gaussian random spheres with a surface 
modulation of cr~0.3 in the case of particulates with irregu- 
lar grain shapes. The main characteristics of these DFFs is a 
strong and relatively sharp peak at L~0. 1 . 

2. DFFs calculated for aggregates of spheres with a fractal di- 
mension D/=2.4 or D/=1.8 in the case of particulates with 
roundish grain shapes. These DFFs peak at L=0.33-0.45 and 
also have a component at very low L values, which is re- 
quired to fit the experimental band profiles and seems to be 
correlated with the agglomeration state in the particulate. 

Our results suggest that DFF models, such as the ones for GRS 
with cr~0.3 and fractal aggregates with Dy-2.4 or Dy-1.8, 
would reflect morphological properties of real particles signif- 
icantly better than CDE or DHS models, so would be useful in 
the simulation of cosmic-dust spectra. This could lead to better 
predictions of infrared band profiles or at least to a valid estimate 
of possible band variations. A necessary condition is that grain 
sizes that are small compared to the wavelength can be assumed 
to dominate the spectrum. 

We recognized that the averaging of spectra obtained by the 
DFF model for the individual crystallographic orientations of an 
anisotropic material is incorrect and leads in some cases to sig- 
nificantly wrong band profiles. This is a problem of all com- 
parable theories and can only be overcome by exact simulations, 
such as a full DDA treatment, which however lacks the statistical 
approach and is difficult to implement in comprehensive astro- 
nomical spectroscopic simulations. Therefore, currently a safer 
way may be to rely on aerosol-measured spectra. 
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Fig. 4. Comparison of measured (solid lines) and simulated extinction spectra for the roundish (left column) and irregular (right) 
particle shapes. DFF models used in the simulations are given in the legends. The forsterite spectra have been divided into two 



